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The so-called /^-deformed oscillator (or /^-oscillator) introduced by A. Jannussis, though possesses 
rather exotic properties with respect to other better known deformed oscillator models, also has good 
potential for diverse physical applications. In this paper, the corresponding /i-Bose gas model based 
on /i-oscillators is developed. Within this model, the intercepts X^^'{K) and X^^'{K) of two- and 
three-particle momentum correlation functions are calculated with the goal of possible application for 
modeling the non-Bose type behavior of the intercepts of two- and three-pion correlations, observed 
in the experiments on relativistic heavy ion collisions. In derivation of intercepts, a fixed order of 
approximation in the deformation parameter ^ is assumed. For the asymptotics of the intercepts 
X^^'{K) and X^'^'{K), we derive exact analytical formulas. The results for X^'^'{K) are compared with 
experimental data, and with earlier known results drawn using other deformed Bose gas models. 
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1. INTRODUCTION 

For the last two decades, diverse models of deformed 
oscillators have received much attention due to both 
their unusual properties, as compared vifith the stan- 
dard quantum harmonic oscillator, and their great po- 
tential for being applied in description of miscellaneous 
physical systems showing essential nonlinearities. The 
diversity of applications ranges from, say, quantum op- 
tics and Landau problem to high energy quantum par- 
ticle phenomenology, modern quantum field theory and 
D-brane (string) theory [1-9]. Among the best known, 
extensively studied deformed oscillators we encounter 
such as the g-deformed Arik-Coon (AC) one [10] and 
Biedenharn-Macfarlane (BM) one [11, 12], as well as 
the two-parameter p,(3'-deformed oscillator [l.'!]. Besides, 
there exists the q-deformed Tamm-DancofF (TD) oscil- 
lator [14, 1-5], also explored though to much lesser ex- 
tent [IG]. Unlike all these, there is very modest knowl- 
edge concerning the so-called /i-deformed oscillator. In- 
troduced in [17] more than a decade and half ago, the 
/i-oscillator essentially differs from the models already 
mentioned and exhibits rather unusual properties [Ls]. 

An important direction of applying deformed oscilla- 
tors (deformed bosons) is the elaboration on their base 
of the respective deformed versions of Bose gas model, 
like the q- or p,q-Bosc gas models in [19-21], [fi, X]. The 
g-deformed extensions of Bose gas are used in modern 
approaches to "'He theory, see e.g. [22]. On the other 
hand, g-Bose and p, g-Bose gas models provide an effi- 
cient ground for successful description of the observed, in 
the experiments on heavy ion collisions, unusual behav- 
ior of two-pion and three-pion momentum correlations. 
As shown in [6, 8, 19, 21, 23] the intercept X^^\K), see 
formula (10) below, of the two-particle momentum corre- 
lation function derived in the deformed Bose gas models 
based on the AC- or BM-type q-deformed oscillators, as 



well as on the p,g-oscillators, lead to a successful model- 
ing of the nontrivial, non-Bose type (such that basically 
differ from bosonic case) properties of pion correlation 
intercepts in the RHIC and CERN SPS experiments on 
relativistic collisions of nuclei. 

Inspired by those results, in this paper we undertake an 
analogous study basing on /i-oscillators ([17], [18]). Cen- 
tral place in our treatment is given to calculation, within 
the ^-Bose gas model involving /x-bosons (/i-oscillators), 
of the two- and three-particle correlation function inter- 
cepts A*-^-* (K) and A*-'^-' (K). For these we establish asymp- 
totically constant type of behavior (i.e., asymptotical ap- 
proaching to a constant given solely by fi). In addition, 
we make some comparison of the obtained results for in- 
tercepts with the existing data. 

The paper is organized as follows. In Section 2 main 
facts about the /(-oscillator are given. For our goals 
we use the concept of the deformation structure func- 
tion [24-20] denoted as (p{N). Here we also sketch some 
unusual properties of the /i-oscillator such as the possi- 
bility of accidental degeneracy of energy levels and the 
non-Fibonacci nature of its energy spectrum. Section 3 
is devoted to calculation of the intercepts of two- and 
three-particle momentum correlation functions. Namely, 
within the /i-Bose gas model based on /i-oscillators we 
calculate these intercepts, for fixed order of approxima- 
tion. Derivation of their large momentum asymptotics 
which is one of our main results is given in Section 4. 
Next 5th Section is devoted to a comparison with some 
experimental data, and with analogous results obtained 
earlier within other (g- or p, q-) deformed Bose gas mod- 
els. The paper ends with concluding remarks. Appendix 
A contains the data for the mean values or thermal av- 
erages of diverse orders (for powers of the number oper- 
ator) , as well as some interesting properties of the coef- 
ficients appearing in the mean values. 



2. SETUP OF THE ^-OSCILLATOR 

In the theory of deformed oscillators it is both con- 
venient and cfScient to use the concept of deformation 
structure function, either in the (p{N) version [24, 25] or 
in the f{N) version [20]. In our present treatment we 
prefer to use the former variant f{N). Each specialized 
structure function determines the particular deformed os- 
cillator model through its corresponding oscillator alge- 
bra, a unital algebra given in terms of the elements a, a^, 
N. These generating elements of the algebra of a par- 
ticular model of linear (usual) or nonlinear (deformed) 
quantum oscillator obey the defining relations 



F{N), 



[N, 



,ti_„t 



[N, a] = -a. 



(1) 



(2) 



From (2) it follows that [iV, aat] = [A^, ata]=0. The rela- 
tion (1) involves F{N), function of the number operator 
TV. For the usual quantum oscillator, a'^a = N, aa^ = N+1 
and F{N) = 1. For general deformed oscillator, like in 
[24, 25], we introduce the deformation structure function 
(p{N) according to the relations 



a^a^ip{N), aa'^ = ip{N + 1). 



(3) 



Due to (1) and (3), the basic commutation relation takes 
the form 



aa^ - a^a ^ F{N) ^ ip{N + 1) - ip{N). 



(4) 



That is, the defining commutation relation is unambigu- 
ously given by structure function. The Hamiltonian H 
will be taken in the standard form H = ^{aa^ + a^a). 
Here and below we put Uoj = I. In terms of structure 
function the Hamiltonian is 



H=^(^<fiN + l) + ^{N)). 



(5) 



We use the g-analog of Fock space. Within it, from the 
ground state obeying 

a|0)=0, A^|0)=0, (OjO) = 1, 

all the n-particlc excited states are generated such that 

N\n)=n\n), ip{N)\n) = ip{n)\n). 

The n-particle normalized state is given by the formula 

(at)" 



\n) : 



V^Pinj 



:|0), (^(7i)! = (p(n)v5(n-l)...v?(l), ^(0)! = 1. 



The creation and annihilation operators a' , a when act- 
ing on the n-particle state give 



1 \n) = \/f(n + l)\n + 1), a\n) = ^/^p{n)\n — 1). 



In what follows we focus on the so called /^-deformed 
oscillator [17] for which the generating elements a, a^ and 
N of the corresponding algebra obey (2) and (4), and the 
defining structure function is 



ipiN)=ip,iN) 



N 



1 + I-lN 



(6) 



where /i is the deformation parameter, /i > 0. With 
this, the basic relation and the energy spectrum for the 
Hamiltonian (5) of /i-oscillator take the form 



aa1-a1a^tPf,{N+l)-ip^{N)-- 



N + l 



N 



l+fi{N+l) 1+fiN 



E„ 



2 V 1 + fi{n + 1) ' 1 + fin 



, (7) 



(8) 



Setting ij. — recovers the algebra and the energies of 
the usual harmonic oscillator. 

Note that the linear structure function if{N) = N, 
which determines the standard quantum oscillator, stems 
as particular case in the limit of deformation parame- 
ter(s), such as q ^ 1 for (/-oscillators or yu — > in the 
present case of ^-oscillator, from each model of deformed 
oscillators (a kind of consistency). 

Unusual properties of fi-oscillator 

The ground state energy and the large n asymptotics 
for the energy of ^-oscillator are respectively Eq = ^,^\_ . 

and Eoo = - . Both of these values basically differ from 
the analogous value and limit for the usual oscillator, 
Eq = ^ and E^ — > oo. Moreover, the ground state 
Eq of the ^-oscillator differs also from that of the Fi- 
bonacci oscillators such as p,g-deformed family and its 
one-parameter q-deformed subfamilies (AC, BM, TD), 
for which zeroth level energy is So = | exactly as for the 
standard oscillator. 

The just mentioned unusual feature of the /i-oscillator 
is tightly connected with the next its peculiarity: the fi- 
oscillator unlike the p,q-deformed one does not belong to 
the class of Fibonacci oscillators. Namely, as shown in de- 
tail in the recent paper [18] the /x-oscillator is typical rep- 
resentative of the wide class of so-called quasi-Fibonacci 
oscillators. 

It is important that the /i-oscillator admits nontrivial 
extensions, in particular through combining it with other 
{q- or p,q-) deformed oscillators. Any such extension will 
inherit the peculiarities of those simpler deformed oscil- 
lators from which it is composed, as shown in [IS]. As 
can be demonstrated for positive fi such that /i < 1, the 
/i-oscillator admits occurrence of accidental energy level 
degeneracy if it is extended by at least one deformation 
parameter. However, unlike the case of Tamm-Dancoff 
q-deformed oscillator [!()], here not all types of pairwise 
energy level degeneracies are possible. 



In other words, direct two- or three-parameter exten- 
sions of //-oscillator can exhibit the phenomena of ac- 
cidental degeneracy (arising in the absence of any ob- 
vious underlying symmetry) of energy levels. In more 
detail this property will be studied elsewhere (see how- 
ever [16, 27-29] where the ability to display accidental 
energy levels degeneracy has been demonstrated for the 
p,g-dcformcd and some g-deformed oscillators). 

3. THE USE OF /i-OSCILLATORS IN /i-BOSE 
GAS MODEL 

Any real system of particles is far from being ideal. 
Therefore, a necessity arises to make complete account 
of inter-particle interactions. Besides, the nonzero proper 
volume or substructure of particles may also be of impor- 
tance. For all these issues it is known that the system 
admits effective description by means of deformed oscil- 
lators [30-32]. Due to that, the system can be treated 
as an ideal one, but consisting of unconventional appro- 
priately deformed "bosons". As result, account of any of 
the indicated complications or reasons turns out to be ef- 
fectively incorporated just in the particular deformation, 
with its specific deformation parameter(s). 

Such an effective picture which uses the g-deformed os- 
cillators has been applied to develop the related g-Bose 
gas model, with subsequent successful application to su- 
perfluid ''He [22] or to an effective description [19-21, 23] 
of the observed non-Bose type behavior of the intercept 
A(m, K) of two-particle momentum correlation function 
of pions produced and registered in the experiments on 
heavy- ion collisions. "Non-Bose" means that the existing 
data for the intercept clearly demonstrate the trend and 
the values that principally differ (see e.g. [33-35]) from 
the would-be constant intercept equal to 1 at any momen- 
tum, for the pions viewed as pure bosons. More advanced 
picture which involves the two-parameter p,(7-deformed 
oscillators and treats n-particle correlation function in- 
tercepts within the p,g-Bose gas model, along with anal- 
ogous application, has been developed in [(>, 8]. Below 
we exploit similar ideas for the case of ^-Bose gas model. 

3.1. Intercept of two-particle correlation function 

In the experiments on heavy ion-collisions, the collid- 
ing (primary) particles generate, as the products of colli- 
sions, a miriads of secondary particles, e.g. pions, which 
are then detected. As known [36-38], correlation func- 
tions carry important information concerning geometry 
and dynamics of the emitting sources. The probability 
for a joint observation of two particles with momenta ki 
and k2 is given by P2{ki,k2). The two-particle momen- 
tum correlation function is defined, see e.g. [36-38], as 



C(2)(fci,fc2)-7 



P2(fcl,fc2) 
Pl{kl)Pl{k2) 



(9) 



izing factor. We put 7 = 1 in what follows. It is 
common to denote the difference of 4-momenta of the 
pair hy Q = ki ~ k2 and the mean 4-momcntum by 
K = {ki + k2)/2. In the case of two identical (of mass 
m) bosons and coinciding momenta, the Q is vanishing. 
Then the two-particle correlation function in variables 
K, Q, and the reduced case Q = of it, is 

C(2)(Q,i^) '"'-^^ C(2)(Q = 0,i^) = 1 + A(2)(m,K) (10) 

where A^^-* (m, K) does appear. The entity A*-^' {m, K) 
at an early stage was usually called the "coherence" or 
(non-) chaoticity parameter. On the other hand, since 
C^^^IQ, K) at fci = fc2 attains its maximal value, X'^'^\K) 
is also termed the intercept of the correlation function 
C'''^\Q, K) (and similarly for 71-th order correlation func- 
tions C*-"^ ) . We will prefer to use the latter term because, 
as stressed in [8] , even in case of complete chaoticity of 
sources the intercept can attain the values differing from 
the would-be purely bosonic ones [39]. The intercept can 
be expressed [3(i] in terms of particle (whose momentum 
is K) creation and annihilation operators as 



A(2)(A') 



a^ a^aa) 
(ata)2 



1, 



K = ki=k2 



(11) 



A natural way to obtain after calculations the nontrivial 
(non-constant) form of the intercept is to use some ver- 
sion of deformed oscillator and to develop on this base 
the corresponding deformed Bose gas model. 

As mentioned, the explicit form of X'^'^^K) has been 
earlier obtained in [6, 19-21, 28] in the framework of de- 
formed versions of Bose gas model based on the popular 
AC-, BM-, p,g-type deformed oscillator models. In the 
case of p,g-bosons this result was extended [(», S] to the 
n-particle correlation function intercepts A'-"''(i^). 

Here we use the set of ^-oscillators to develop the cor- 
responding //-Bose gas model. Within it, we derive the 
intercepts of the two-particle and three-particle correla- 
tion functions. For that goal we take (11) as our starting 
point and in view of (3) make account of the defining 
structure function of the chosen /i-oscillator, namely 



ita = ip^{N) = [N]^ = - 



N 



fiN 



(12) 



where Pi(fci), i = 1,2, means pure single-particle prob- 
ability distribution for i-th particle and 7 is normal- 



From eq. (2), for a function f{N) the property follows: 
af{N) = f{N+l)a, a^ f (N) = f {N - l)al (13) 
From (12) wc have one- and two-particle distributions 

(atfltaa) = {a^N]^a) = {a^a[N-l]^) = {[N]^[N^1]^). (14) 

Let us point out the principal difference between de- 
formed and usual, non-deformed, approaches. Due to 



non-linear relation (12) of a^a and the operator N (for 
the usual oscillator, we have a'^a = N), one-particle dis- 
tribution is given by (a^a), and not by (N). 

Account of the both relations (14) in eq. (11) implies 



a'^HaO 






1. 



(15) 



Assuming that /i is positive and sufficiently small we 
present the ratio (12) as 

oo 

[N]f, = N - fiN^ + f?N^.. =Nj2i~lJ'Ny ■ (16) 



Er. 



-,— /3cjn 



Z-/n 



^~/3ujn 



al3uj _ I 



(19) 



where f3 = -^, k — 1. Since in what follows we will 
use the developed model for description of the intercepts 
of pion correlation functions, we assume the particle's 
energy to be w = [m^ + K^)-'-/^. Similarly to (19) we 
can derive the mean values {N''), k > 2. Let us present 
those of them which are needed in (18): 

/,r3N 1 6 6 



{N)^-, {N^) = -+ ^ 



2 ' 



... 1 14 36 24 

(^ ) = - + — + — + —' 



After substitution of the expression (16) in eq. (15) the 
formula for intercept of two-particle momentum correla- 
tion function reads: 

(17) 
Let us note that we cannot obtain single Taylor series 

for the ratio A/i (/-f) in (15), since the numerator and 
denominator imply their own averages and, besides, we 
have no general formula for necessary coefficients. 

So, in calculating A/i ' {K) we deal with series and hence 
use some approximation. Let us restrict ourselves with 
the approximation in (16) up to ^^ . In this case the last 
formula (17) can be cast to the form 



a('Ha')U5 = 



ELia.(^'^) 



(ELo(-i)^m^^-'+^)' 



-1 



(18) 



with the coefficients a^ as follows: 
Qfi = — 1 — /i — /i^ — /i'^ — ^"^ — ^^ , 

"2 = 1 + S^t + 4/i2 + 5/1^ + 6/i'' + 7^^ 
as = -2^ - 6m^ - lO/i^ - I5;i4 _ 21^5, 
UA = 3/i2 + lO/i^ + 20/i^ + 35/i^ 
as = -4^3 _ it^^i _ 35^5^ 

ag = 5/i** + 21//^, 
ar — — 6/i'''. 

Next step of obtaining Aj^ (K) consists in finding mean 
values of powers of TV. Say, (N) yields (like in [40]): 



(N) 
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The expressions for mean values of these and higher order 
powers TV™ are of general form 



(iV")==^B(" 



ofSi^ 



1. 



(20) 



The coefficients Br , both for tti < 7 and for 8 < m < 
14, see Appendix A, show very interesting properties, 
also mentioned in Appendix A: say, they are arranged in 
a remarkable analog of the well-known Pascal's triangle. 
To summarize: the explicit expression for A^^' {m, K) 
in the /i^ approximation is given by eq. (18), with the 
above mean values {N^). This result covers the lower n^ 
and /i^ orders as well. Using data in Appendix A, similar 
result can be given for the order /i^, and so on. 

3.2. Intercept of three-particle correlation function 

Like in the preceding Subsection, the intercept of three- 
particle momentum correlation function is expressible in 
terms of creation and annihilation operators: 

^i^K) ^ /„t„v.s - 1' A' = fci =k2^ ks. (21) 



(ata) 



fa)3 



(3), 



Using (13) and (12) the intercept A}i (K) takes the form 

^(3).^^ {[NUN-lUN-21) 

Account of (16) in (22) leads to general formula for in- 
tercept of three-particle momentum correlation function: 



,(3).,.^ _ {NiN~l){N^2)E^=o{~^^NrEn=oi-^^iN-l)rEZoi-^^iN-2)Y) 



(23) 



Like before, we consider the order /z^ approximation and 
obtain the expression for A/^ {K) in the form 



A(^)WL= 



-1 



(E'=o(-i)>'^^'^+') 

where the coefficients r\k are the foUowing: 
?7i = 2 + 6pt + 14/i2 + 30/1^ + 62/1^ + 126^^ 
7^2 = -3 - 15/i - 45/i2 - 115m^ - 273Ai* - 623^'5, 
?73 = +1 + 12/i + 55/i2 + 180/i3 + Sll^i** + 1344^^ 
7^4 == -3/i - 30/i2 - 145^^ - 525^^"* - 1652/x^ 
% = 6/i2 + 60// + 315^"* + 1260/l^ 
% = -lO^i^ - 105Ai* - 602/i^ 
r/7 = 15m'' + 168//^ 



(24) 



'78 



-21^5^ 



4. ASYMPTOTICS OF INTERCEPTS OF TWO- 

AND THREE-PARTICLE CORRELATION 

FUNCTIONS 

Here we study the asymptotical behavior at /3w — ?> oo 
(i.e. at |K| -^ oo, constant T) of the intercepts A/^ and 
A/i of the two- and three-particle correlation functions. 

Let us find the asymptotical expression for A/^ within 
the order [i' approximation, see (18). With the account 
of relevant averages from Appendix A both in the nu- 
merator and in the denominator of (18), we derive the 
result 



Ai'2(M)l, 



1 - 4/i + \\[P- - 26jj^ + 57//-* - 120/i^ 



1 



2/i + 3/i2 - 4/^3 + 5Ai4 - 6/^5 
This expression can be written in the form 



(l + A.)2(l + 2/0V(l + Ai)2 



1 + 2^ 



(27) 



It should be noted that in the asymptotics the inter- 
cept does not depend on the temperature, nor on the 
momentum and mass of particle. It is the deformation 
parameter /x alone that determines the intercept asymp- 
totics. So Xas. (a*) is a number < 1 at any fixed /i > 0. 

Now consider the asymptotics of the three-particle cor- 
relation function intercept. Applying the limit /?aj — >■ oo 
to eq. (24), with the {N'^) taken from Appendix A, we 
obtain the desired asymptotic formula in fi^ approxima- 
tion: 



6(1 -6/i + 25/^2 - 90/i3 + 30lAJ**-966Ai5) 



(28) 

Using this expression (which also contains the order fi^, 
/z* results) we find the exact asymptotical formula 



aS(m) 



6Er=o(-ir(5(3^+'+i)-2^+' 



Er=o(-i)i(^+i)(^+2)M'^ 



(1+M)(1+2/.)(1+3m) 



- 1 = 



7^^ 



(T+mF 



(l + 2/i)(l + 3/i) 



(29) 



In fig. 1 we present the two- and three-particle correlation 
function intercepts Ajj , Aji for the order /i^ approxima- 
tion, with their asymptotes shown explicitly. 
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FIG. 1. Left: The intercept A'^^ to order /i^, with asymptotes 

(2) 

(dashed lines). Asymptote for /.i=0.1 curves is XaJ. ~ 0.8331 
and for fi—0.l5 curves is AL = 0.7664. Right: The intercept 
A^ ' to the order f/" , with asymptotes. Asymptote for /i=0.1 

2.9964. 
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curves is XaJ. = 3.6365 and for fi—0.l5 curves is AJi 



(3) 



5. INTERCEPTS IN /i-BOSE GAS MODEL 

COMPARED WITH OTHER MODELS AND 

EXPERIMENTAL DATA 

It is of interest to compare the behavior of the inter- 
cepts A^^^ and A'^' obtained in the /x-Bose gas, see fig. 1, 
with main features of such intercepts found within g-Bose 
gas models [8, 19-21, 28] and plotted in the corresponding 
figures of those works. As seen from the comparison, the 
picture is qualitatively the same: with increasing mean 
momentum, in each model we observe the occurrence 
of minimum, convexity downwards, then the transition 
through inflection point to the convexity upwards, and 
finally asymptotical tending to a constant value strictly 
lesser than 1 if g 7^ 1 or /i 7^ 0. 

Now consider some experimental data, e.g., taken from 
[33]. It is nice that the experimental data for the two- 
pion correlation intercepts (4 points corresponding to 4 
momentum bins) respect the basic features mentioned in 
the previous paragraph (the data show stable qualitative 
picture or "trend", for different colliding nuclei, energies 
of collision and centralities) . In our fig. 2 (left), for il- 
lustration, we make a comparison of the curve A^^^ from 
(18) with experimental points [■'>'■-'>] corresponding to the 
most central 62.4 GeV Au+Au collision. In principal, 
we could get besides qualitative, also better quantita- 
tive agreement with data within the considered /i-Bose 
gas model, although by taking lesser temperature and 
slightly greater value of fi. Note again that the fig. 2 
{left) serves for illustration only, and more detailed and 
realistic comparative analysis will be given elsewhere. 

What concerns the existing data for three-pion corre- 
lation intercept, these are obviously insufficient. For that 
reason we postpone confronting with data of our analyt- 
ical results on A^'^-' (A') for future work as well. 

Before ending this section, let us note that an impor- 
tant function r^^'{K) was introduced in [41]. It is com- 
posed from A*^^' (K) and A'^' (K) according to the formula 



,(3).j.^ l A(3)(A-)-3A(^)(i^) 

^ ' 2 (A(2)(/v))3/2 ■ 



(30) 



Due to special form of the ratio, the function (30) pro- 
vides improved purity (on the experimental side). In- 
deed, many unwanted distortions which especially for low 
K may affect \^^'>{K) and A'^'^^(-ftr) when considered sep- 
arately, in r^^) {K) are mutually canceled. 

In case of g-Bose gas this function was studied in [8] 
by one of us. From our analytical results (18) and (24) 
derived in the /x-Bose gas model, the formula for the func- 
tion rjj (K) reads 



y:L,vi{'^') 



1:1=, ^^N") 



r(3) 



i'>iK) 



(eLo(-i)>=^'=+0' (i:"=o(-i)>'-jv'-+i) 



i:Li"'^<^^''> 



(e"=o(-i)>'-^'^+0 



T-i 



3/2 



Finally, using exact asymptotical formulas (27) and (29) 

(3) 
we find the exact expression for the r), (K) asymptotics: 



„(3) 



(m) 



1 + 3^1 



VT 



i^i. 



(32) 



Fig. 2, right, gives the plot of r\i'{K) along with its 
asymptotics. As seen, small raise of the parameter /i 
results in basic change of the properties of r), {K): from 
positive to negative sign. But, as seen from (32) its 
asymptotical value is always positive since /i < 1. 
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FIG. 2. Left: The intercept A{f'(A') from (18), to /i^ versus 
experimental data [•")•)]. The parameter values are ^ = 0.255, 
T = 100 MeV. Right: The function r'}^\K) from (31), to /x^ 
with its asymptotes (dashed horizontal lines). Asymptote for 



.(3) 



the /i=0.1 curves is raJ. = 0.7472 and for the /i=0.15 curves 



(31) 



is rS = 0.5194. 



6. CONCLUSIONS 

The /i-deformed oscillator model is a peculiar one pos- 
sessing novel interesting features with regards to the 
other, more popular and better known models of de- 
formed oscillators. The specifics of the /i-oscillator is 
rooted in its defining structure function whose non- 
polynomial, rational dependence on the number operator 
A'' causes nontrivial properties of this oscillator especially 
those mentioned at the end of sect. 2. 

The energy spectrum of /i-oscillator docs not form [18] 
the usual Fibonacci sequence. For that reason we treated 
the /i-oscillator as the one obeying a generalization of 
Fibonacci relation and called it quasi-Fibonacci oscilla- 
tor [1^]. We hope that the deformed oscillators possess- 
ing quasi-Fibonacci property may be useful not only in 
the present context of multiparticle correlations of pions 
(kaons) produced and registered in the experiments at 
RHIC or LHC, but also, say, for description of the prop- 
erties of some quasi-periodic chains [42] . 

In the present paper, the /i-oscillators are employed for 
elaborating on their base the corresponding /i-Bose gas 
model. The results concerning one-, two-particle distri- 
bution functions and the two-particle correlation func- 
tion intercept obtained in /t-Bose gas model are derived 
by using definite (/i^) order of approximation. From the 



obtained intercept of two-particle correlation function de- 
pending on the mean momentum K, the exact analytical 
formula for its asymptotics is found. We observe the im- 
portant fact that in the asymptotics, the intercept does 
not depend on the particle mass, nor on the tempera- 
ture or momenta: the deformation parameter ^ is the 
only entity that survives and completely determines the 
/3w — !> oo asymptotics of A'^-'(/i). Analogous set of re- 
sults, i.e. the explicit formulas and the exact asymptot- 
ical expression have been obtained for the intercept A*^^' 
of three-particle correlation function. 

Let us emphasize ones more that the obtained analyt- 
ical formulas for \i,'{K), Xi,'{K) and ri,'{K) provide 
explicit and nontrivial dependence on the particles' mean 
momentum. This allows to make a comparison with ex- 
perimental data on intercepts of pions, as shown in spe- 
cial case, see fig. 2, left. Certainly, more detailed data are 
needed in order to make choice in favor of a particular 
model and we believe the situation will soon improve. 

We hope that the above /x-Bose gas model constructed 
by employing /x-oscillators will find interesting applica- 
tions in different branches of nonlinear quantum physics: 
in description of high-energy quantum physics phenom- 
ena, in *He theory, in quantum optics, etc. 

ACKNOWLEDGMENTS 

The authors are thankful to the referees for valuable re- 
marks and recommendations. This research was partially 
supported by the Grant 29.1/028 of the State Foundation 
of Fundamental Research of Ukraine and by the Special 
Program of the Division of Physics and Astronomy of the 
NAS of Ukraine. 



[1] V.I. Man'ko, G. Marmo and F. Zaccaria, Phys. Lett. A 

197, 95 (1995). 
[2] G. Alvarez ei al, J. Phys. A: Math. Gen. 35, 8705 (2002). 
[3] J. Crnugelj, M. Martinis and V. Mikuta-Martinis, Phys. 

Rev. A 50, 1785 (1994). 
[4] J. Ben Geloun, J. Govaerts and M.N. Hounkonnou, Eur. 

Phys. Lett. A 80, 30001 (2007). 
[5] A. Algin, J. Stat. Mech.: Theor. Exp. 04007 (2009). 
[6] L.V. Adamska and A.M. Gavrilik, J. Phys. A: Math. Gen. 

37, 4787 (2004). 
[7] C.L Ribeiro-Silva, E.M.F. Curado and M.A. Rego- 

Monteiro, J. Phys. A: Math. Theor. 41, 145404 (2008). 
[8] A. M. Gavrilik, SIGMA 2, Paper 074, 12 pages (2006). 
[9] K. Okuyama, Phys. Lett. B 673, 232 (2009). 
[10] M. Arik and D.D. Coon, J. Math. Phys. 17, 524 (1976). 
[11] L.C. Biedenharn, J. Phys. A: Math. Gen. 22, 873 (1989). 
[12] A.J. Macfarlane, J. Phys.A: Gen. 22, 4581 (1989). 
[13] A. Chakrabarti and R. Jagannathan, J. Phys.A: Math. 
Gen. 24, 711 (1991). 



[14 

[15; 

[16 

[17 
[18 

[19 

[20 

[21 

[22; 

[23; 

[24; 

[25; 

[26; 

[27 

[28; 
[29; 

[30 

[31 

[32' 
[33; 

[34 
[35| 

[36 

[37 

[38; 
[39 

[40 



[41 
[42; 

[43] 



K. Odaka, T. Kishi and S. Kamefuchi, J. Phys. A: Math. 
Gen. 24, 591 (1991). 

S. Chaturvedi, V. Srinivasan and R. Jagannathan, Mod. 
Phys. Lett. A 8, 3727 (1993). 

A.M. Gavrilik and A.P. Rebesh, Mod. Phys. Lett. A 22, 
949 (2007). 

A. Jannussis, J. Phys. A: Math. Gen. 26, 233 (1993). 
A.M. Gavrilik, I.I Kachurik and A.P. Rebesh J. Phys. A; 
Math. Theor. 43, 245204 (16pp) (2010). 
D. Anchishkin, A. Gavrilik and N. lorgov, Eur. Phys. J. 
A7, 229 (2000). 

A.M. Gavrilik, Nucl. Phys. B Proc. Suppl. 102, 298 
(2001). 

D. Anchishkin, A. Gavrilik and N. lorgov, Mod. Phys. 
Lett. A 15, 1637 (2000). 

M. Rego-Monteiro, L.M.C.S. Rodrigues and S. Wulck, 
Physica A 259, 245 (1998). 

Q.H. Zhang and S.S. Padula Phys. Rev. C 69, 024907 
(2004). 

S. Meljanac, M. Milekovic and S. Pallua,Phys. Lett. B 
328, 55 (1994). 

D. Bonatsos and C. Daskaloyannis, Prog. Part. Nucl. 
Phys. 43, 537 (1999). 

V.i. Man'ko, G. Marmo, E.C.G. Sudarshan and F. Zac- 
caria, Phys. Scripta 55, 528 (1997). 

A.M. Gavrilik and A.P. Rebesh, Ukr. J. Phys. 53, 586 
(2008). 

A.M. Gavrihk and A.P. Rebesh, Mod. Phys. Let. A 23, 
921 (2008). 

A.M. Gavrihk and A.P. Rebesh, SIGMA 3, Paper 112, 
11 pages (2007). 

A.M. Scarfone and P. Narayana Swamy, J. Phys. A: 
Math. Gen. 41, 275211 (2008). 

S.S. Avancini and G Krein, J. Phys. A: Math. Gen. 28, 
685 (1995). 

W.R. Perkins, Int. J. Th. Phys. 41, 833 (2002). 
B.I. Abelev et al. (STAR CoUab.), Phys. Rev. C 80, 
024905 (2009). 

T. Humanic, Int. J. Mod. Phys. E 15, 197 (2006). 
D.V. Anchishkin, A.M. Gavrilik and S.Y. Panitkin, Ukr. 
Phys. J. 49, 935 (2004). 

M. Lisa, S. Pratt, R. Soltz and U. Wiedemann, Ann. Rev. 
Nucl. Part. Sci. 55, 357 (2005). 

U.A. Wiedemann and U. Heinz, Phys. Rep. 319, 145 
(1999). 

T. Csorgo, Heavy Ion Physics 15, 180 (2002). 
Recall that for ordinary bosons the intercept A"' should 
be constant: A^^^^ — Cg^^^ — 1 = 1. Note also that m in 
(10) is the mass of (identical) registered particles. 
N.N. Bogoliubov and N.N. Bogoliubov, Jr., Introduc- 
tion to Quantum Statistical Mecfianics (Nauka, Moscow, 
1984), 384 pp. 

U. Heinz and Q.H. Zhang, Phys. Rev. C 56, 426 (1997). 
Luca Dal Negro and Ning-Ning Feng, Optics Express 15, 
14396 (2007). 

N.J. A. Sloane and S. Plouffe, The encyclopedia of integer 
sequences, Academic Press (1995). 



J 



APPENDIX A. EXPLICIT EXPRESSIONS FOR THE MEAN VALUES (iV™) 



(m) 



.;3cu 



The expressions for mean values of the powers iV™ are given as (N™) — J2T=i ^ 
coefRcients Br , m — 8, ..., 14, placed in the Table 1. For m = l, ..., 7 see the expressions preceding eq. (20). 

TABLE 1. The coefRcients Br' for mean values of the powers iV™ (8 < m < 14) 



1, with the 



r\m 


8 


9 


10 


11 


12 


13 


14 


1 


1 


1 


1 


1 


1 


1 


1 


2 


254 


510 


1022 


2046 


40946 


8190 


16382 


3 


5796 


18150 


55980 


171006 


519156 


1569750 


4733820 


4 


40824 


186480 


818520 


3498000 


14676024 


60780720 


249401880 


5 


126000 


834120 


5103000 


29607600 


165528000 


901020120 


4809004200 


6 


191520 


1905120 


16435440 


129230640 


953029440 


6711344640 


45674188560 


7 


141120 


2328480 


29635200 


322494480 


3162075840 


28805736960 


248619571200 


8 


40320 


1451520 


30240000 


479001600 


6411968640 


76592355840 


843184742400 


9 





362880 


16329600 


419126400 


8083152000 


130456085760 


1863435974400 


10 





— 


3628800 


199584000 


6187104000 


142702560000 


2731586457600 


11 





— 


— 


39916800 


2634508800 


97037740800 


2637143308800 


12 





— 


— 


— 


479001600 


37362124800 


1612798387200 


13 





— 


— 


— 


— 


6227020800 


566658892800 


H 





— 


— 


— 


— 


— 


87178291200 



The coefficients Br in (-^*) o,nd an analog of Pascal's triangle 

Here wc examine in some detail the coefficients appearing in the mean values of powers of N, which determine the 
intercepts A^^^ or A^"^^. Let us specialize (20) for to = 1, 2, ..., 7. Then, e.g., for m = 5 we have: B\ ' = 1,B2 - 



?(5 



:30,Bf^: 



(5) 



150, i?4"''' =240,i35"''' = 120. Consider the correspondence between the coefficients Bi- in (A^™) and the entries of 
some analog of Pascal's triangle. For the classical Pascal's triangle [l-S], there is a "rule": each coefficient in m-th row 
equals the sum of two adjacent coefficients in the (to — l)-th row. Now arrange our coefficients as the entries of some 



(m) 



1, ...,m, for seven rows (m < 7) looks as 



triangle. Then, the triangle with entries B,.' 

1 
1 2 

1 6 6 

1 14 36 24 

1 30 150 240 120 

1 62 540 1560 1800 720 
1 126 1806 8400 16800 15120 5040 

With regards to the famous Pascal's triangle, here we have two essential distinctions. First, now we encounter unities 
only on the left hand side. Second, the above mentioned "rule" is to be modified. That is, each entry Br in 77i-th 



row equals the sum of two adjacent entries above it, multiplied by r, the position number of Br 



('"). 



b: 



(m) 



r fi 



(m-l) 



B. 



(m-l)^ 



It is worth to note that in the mean value (iV™), the last coefficient fixed by ?' = ?Ti is given by the formula 

Btl = ml . (33) 

Here to is both the power of N and the row number, while r is the number of entry in m-th row. Likewise, the entries 
with r = 7Ti — 1 in the m-th row, and the second entry of each m-th row are given as 



B, 



(m) 



(to-1) 



B, 



(m) 



= 2(2" 



- 1) = 2" - 2. 



(34) 



Let us note that (33) and the first expression in (34) coincide with the Lah numbers L{m^ 1) and i(m, 2) respectively, 
while the second formula in (34) coincides with the Mersenne numbers (see [43]). 
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